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ABSTRACT 

Approximation formulas are found forx(O) andx(l), where 
x(t) satisfies 


x - f(x, t) + A(t)x , 
x(0) = x Q1 x(l) = x x . 
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APPROXIMATE SOLUTION OF A TWO-POINT 
BOUNDARY VALUE PROBLEM 


INTRODUCTION 

Consider the following boundary- value problem: 


x - f(x, t) + A(t)x , 
x(0) = x 0 , x(l) = x, , 


( 1 ) 


where t is a scalar, x is a column matrix, A(t) is a matrix function of t, and 
the overdots indicate differentiation with respect to t . 

Formulas will be found for x Q and x x such that 

x 0 = 0(0), i x = 0(1) , 


where 0(t ) satisfies 

0(0) = x(0), 0(0) = x(0), 0( 0 ) = x(0) , 

0(1) = x(l), 0(1) = x(l), 0(1) = x(l) . 

It will not be necessary to find 0( t ) but only to make certain assumptions as to 
its form. 

The author has found the method useful in the calculation of transfer tra- 
jectories for space vehicles and in preliminary orbit determination from 
observations of the position of a spacecraft at two times. 


AN APPROXIMATE SOLUTION 

Let a be any component of x, a(0) =a 0 ,a(l) = a x , and assume approxi- 
mations of the form 



(3) 


<*1 - c o a 0 + d O a l + c 2 a 0 + d 2<*l + C 3 s 0 + d 3 S l • 

The scalars a Q , b 0 , a 2 , b 2 , a 3 , and b 3 are determined by assuming (2) to be 
exact when a = <p i (t ) , i = 0, • • *, 5, where the cp^s are linearly independent 
over the interval [o, l] with derivatives through the third order at t = 0 and 
t = 1 . The coefficients in (3) are determined in a similar way. If we use the 
same set of </>.'s for each component of x, we can write 


*0 " a 0 X 0 + b 0 X l + a 2 X 0 + b 2 X l + a 3 X 0 + b 3 x i > 


(4) 


x, - c 


0 X 0 + d 0 x i + C 2 X 0 + d 2 X l + C 3 X 0 + d 3 X 1 • 


(5) 


We eliminate the third derivatives in (4) and (5) by using 


x - (P + A)x + Ax + w , 


( 6 ) 


obtained from (1), where P is a matrix with element in the ith row and j th 
column equal to the value of 9f 1 /dx* t and w is a column matrix with ith ele- 
ment equal to the value of d f Vat, f 1 and X J being respectively the ith 
component of f and the jth component of x . Substituting (6) into (4) and (5), 
we obtain 


i 1 ' a 3 Q 0 ) x o " b 3 Qi x i “ z 3 - 

(7) 

" c 3 Qo x o + i 1 ~ d 3 Q i) x i = y • 

(8) 


where I is the unit matrix and 


Q = P + A , 


y8 

y 


a o x o 

+ 

Vi 

+ 

K 1 

c o x o 

+ 

d 0 x l 

+ 

( C 2 J 


+ a 3 A o ) X 0 + ( b 2 J + b 3 A l)*l + a 3 W 0 + b 3 W l • 
+ C 3 A o) X 0 + { d 2 1 + d 3 A l) X l + C 3 W 0 + d 3 W l * 
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Solving (7) and (8) we obtain 


(B + bQ^i, 

= P + Q,(v - d 3 ^)> 

(9) 

(B + bQ 0 Q,)i, 

= y + Q 0 ( c 3^ ~ a 3 y)- 

(10) 


B I “ a 3^0 ~ d 3^1 ’ 

b ~ a 3 d 3 ” ^3 C 3 ' 


A POLYNOMIAL APPROXIMATION 
Substituting successively 

a - <p i ( t) = t l , i “ 0, 1, 2, 3, 4, 5, 
into equation (2), we obtain 


b Q + 2a 2 + 2b 2 = 0 , 

b Q + 6b 2 + 6a 3 + 6b 3 = 0 , 

b 0 + 12b 2 + 24b 3 = 0 , 

b 0 + 20b 2 + 60b 3 = 0 . 
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The solution of this set of equations is 


= -1, b n - 1, a. 


7 

20 


3 

20 


" 20 ’ b 2 


In a similar way we find 


c n ”1 > dft ^ > c ' 


_L . 

20 ’ d 2 


7 

20 


1 

30 


1 

20 


30 • 
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